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We study a mass transport model, where spherical particles diffusing on a ring can stochastically 
exchange volume v, with the constraint of a fixed total volume V — ^ being the total 

number of particles. The particles, referred to as p-spheres, have a linear size that behaves as v\^^ 
and our model thus represents a gas of polydisperse hard rods with variable diameters v\^^ . We show 
that our model admits a factorized steady state distribution which provides the size distribution 
that minimizes the free energy of a polydisperse hard rod system, under the constraints of fixed A'^ 
and V . Complementary approaches (explicit construction of the steady state distribution on the 
one hand ; density functional theory on the other hand) completely and consistently specify the 
behaviour of the system. A real space condensation transition is shown to take place for p > 1: 
beyond a critical density a macroscopic aggregate is formed and coexists with a critical fluid phase. 
Our work establishes the bridge between stochastic mass transport approaches and the optimal 
polydispersity of hard sphere fluids studied in previous articles. 
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I. INTRODUCTION 



Condensation phenomena in stochastic models of mass transport is a subject of much current interest, for reviews 
see [l|, 0] ■ Typically, in these models a globally conserved quantity, which for the purposes of this introduction we refer 
to as mass, is transferred stochastically between sites of a lattice according to some prescribed dynamical rules. One 
is interested in the properties of the stationary state generated by the stochastic dynamics, for example the single-site 
mass distribution which is the probability distribution of the amount of mass at a lattice site. Such models provide 
both microscopic and effective descriptions of the dynamics of various complex systems, for example traffic flow [1, 
granular clustering phase ordering @, network rewiring 0,@|, force propagation [oj, aggregation and fragmentation 
[l3,[IH and energy transport p^ . In this paper we study another realisation of stochastic mass transport in a new 
context: the sampling of polydispersity in simulations of hard-rod systems. 

Condensation occurs in stochastic mass transport models, when above a critical value of the global density a finite 
fraction of the mass condenses onto a single site d, [13, [13, [l3| • In some ways the phenomenon is similar in character 
to Bose-Einstein Condensation, but in contrast to Bose-Einstein Condensation, it occurs in real space i.e. at a lattice 
site. The signature of condensation is seen in the steady-state single-site mass distribution. Below the critical density 
the distribution typically decays exponentially for large mass, indicating a fluid phase. At the critical density the decay 
of the single-site mass distribution is slower^ typically it decays as a power law or sometimes a stretched exponential 
distribution, indicating a critical fluid jisl . Above the critical density a bump in single-site mass distribution 
emerges and corresponds to a single site containing the excess mass above the critical value. Thus the condensed 
phase consists of a condensate co-existing with the critical fluid. 

In the present work we consider the stochastic dynamics of mass transfer in a rather different context, that of 
sampling polydispersity in a hard rod fluid. Polydispersity arises naturally in a wide variety of natural and synthetic 
materials. In some cases, polydispersity is an intrinsic property of the system; the elementary constituents are 
characterized by a varying physical property (size or charge, for example) which does not change as a result of the 
particles' interactions. A different type of polydispersity arises as a result of the interactions among constituent 
elements in situations where they self-assemble to form aggregates of varying size, for example. In this second 
situation, the flnal size distribution emerges as a result of the interaction among the constituent elements and it will 
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be controlled by the dynamical constraints and basic symmetries underlying the kinetic processes [17|, [l^ . Previously, 
in [l^ a stochastic algorithm has been used which allowed hard spheres to diffuse and to exchange volume, subject to 
the hard core constraint. It turned out that above a certain volume fraction small numbers of large particles would 
emerge in the distribution generated by the dynamics. Within approximate theories of the Percus-Yevick family, this 
was understood as a continuous phase transition ^17l] . 

In this paper, we will show that the transition observed in polydispersc hard spheres can in fact be understood 
as a condensation transition arising from constraints in the configuration space. To make the problem analytically 
solvable, we somewhat simplify matters by constraining the spheres to move on a one-dimensional ring (so that the 
particles effectively behave as hard rods). Basically there are two conserved quantities in play: the total volume of the 
spheres and the linear size (length) of the lattice. However, a new feature which appears, compared to previous studies 
of condensation in systems with two conserved quantities [20| , is a configurational constraint which effectively couples 
the two conserved quantities. The configurational constraint results from the hard-sphere condition and requires that 
X, the separation of a sphere from its neighbour must be greater than the diameter of the sphere defined as £ = v^^p 
where v is the volume of the sphere {xi > £i in the notation of Fig. [T]) . We retain p as a parameter of our particles 
which we refer to generally as p-spheres: for discs p = 2 and for spheres p = 3. The dynamics comprises diffusion, 
which is a stochastic exchange of inter-particle distances, and in addition there is a stochastic exchange of volume. 
When the volume fraction is large the dynamics becomes constrained. For p > 1, beyond a critical volume fraction, 
condensation occurs. That is, in the stationary state one large sphere emerges containing the excess volume and 
leaving the rest of the system in the critical fluid phase. When p < 1 there is no condensation but a vestige of the 
transition remains wherein above some density threshold, the volume distribution develops a bump around a finite 
characteristic volume, unlike the low density behaviour where the most probable volume is 0. 

We develop two complementary approaches to understand the condensation phenomenon: one microscopic and 
one thermodynamic. In the former approach we define a microscopic dynamical model for which we can solve the 
equilibrium state exactly. We then use the usual machinery of statistical mechanics to solve the model. In the latter 
approach we write down a free energy functional, and then use scaling arguments to deduce various thermodynamic 
relations, and arrive at the same form for the equilibrium distributions as with the microscopic approach. 

A. Summary of results 

As a guide to the reader, we first summarise how the paper is organised and identify the main results obtained from 
our calculations. In section |TT] we define our model and determine simple microscopic dynamical rules which lead to 
an equilibrium distribution of the factorised form Q. In particular we identify rules which lead to an equilibrium 
with equal probabilities for all allowed configurations pH)) and in section ITIIl we consider the thermodynamics of this 
equilibrium state. The calculations are carried out within the grand canonical ensemble for the system which we 
define in IIII Al and within which the marginal distributions for volume and separation of a particle p{v) and p(x) 
are determined (|18ll9p . These distributions are expressed in terms of the two Lagrange multipliers corresponding to 
the constraints of mean density p and volume fraction (f> (|23p . In section IIIIBI we introduce a single combination of 
these two Lagrange multipliers which allows us to fix p and consider the reduced volume fraction (f)* , defined in (|35p , 
as a function of a single variable p6|) . We define the entropy (p9)) in section IIII CI and show that the entropy has a 
maximum as a function of (p* . We then analyse the two cases p < 1 and p > 1 separately in sections IIIIDI and IIII El 
showing that for p < 1 the entropy increases to a maximal value then decreases again as (j)* is increased, where for 
p > 1, the entropy sticks to its maximal value when (p* is increased past its critical value given by The latter 

scenario is explained as a condensation transition in IIII El In section IIIIFI we show how the limit of p — > oo recovers 
the results of the Tonks gas. 

In section ITVl we analyse the nature of the condensate by working within the micro-canonical ensemble. By invoking 
results for the large deviations of sums of random variables we show that in the condensed phase the marginal 
distribution for the volume of a particle p{v) has two pieces (|88l89p . The first piece represents a critical fluid 
distribution and the second represents a Gaussian peak corresponding to the condensate. 

Having shown that our model, initially defined through a set of dynamical rules, admits a steady state that is in 
fact the equilibrium state of a hard rod system, we are then in section |V] in position to approach the problem by a 
complementary free energy functional. We show that beginning from the form (|90p for the free energy as a functional 
of the particle size disrtibution, simple considerations imply that the distribution of the particle diameters should 
follow (|100p . We then revisit the condensation transition in section IV CI and show the results of section [III El mav be 
recovered. In section IVTl we present numerical results, based on Monte Carlo simulations of the microscopic model, 
confirming our analytical predictions. In particular. Figure [7] confirms the emergence of a condensate at the critical 
volume fraction, and Figures [51 [5] confirm the predictions (|88l89p of section IIVI for the particle volume distribution. 
We conclude with an overview in section IVIII 
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II. MODEL 



In this work we introduce a simple one-dimensional model of a fluid of hard p-spheres with stochastic dynamics 
comprising diffusion and volume exchange. We show that appropriate choices of the dynamical rates allow one to 
obtain the steady state exactly as a factorised form, in particular the steady state may be such that any allowed 
microscopic configuration appears with equal probability. 

The model consists of p-spheres diffusing on a one dimensional ring and exchanging volume with hard core interac- 
tions (see Figure [Ij . The distance between the left hand side of p-sphere (particle) i and p-sphere z -I- 1 is given by Xi 




FIG. 1: Schematic view of three p-spheres on a line, with definition of the left to left inter-particle distance Xi and gap gi. 



and the diameter of p-sphere i is given hy £i = where Vi is the volume of p-sphere i. The volume of the p-sphere 
with unit diameter has been set to 1 but one could easily generalise to the volume of the unit diameter p-sphere being 
a. 

The model is equivalent to N sites i = I, . . . , N each with site variables Xi, Vi and with periodic boundary conditions 
xjv+i = xi, ujv+i = vi. In the micro-canonical ensemble we have fixed total volume 

Y^v^^V (1) 

i 

and fixed total length L 

Y^x,=L. (2) 

i 

The hard-core interaction between p-spheres implies that for each separation xi we have the constraint 

X. > v\l^ . (3) 

The model then can be interpreted as a system of N polydisperse hard rods on a ring where each rod has a variable 
diameter li = and satisfies the hard rod constraint, namely the gap between successive rods gi > for each i. 
This is a generalization of the classical monodisperse hard rod system, the Tonks gas [2lj , where each rod has the 
same diameter i independent of i. Note that our polydisperse model reduces to the monodisperse case in the large p 
limit since £i = — > 1 (independent of i) in the p — s- cx) limit. 



A. Dynamics 

We consider diffusion and volume dynamics implemented by microscopic transition rates on the variables {xi,Vi}. 
The diffusion dynamics implies dynamics for the separations {xi}. We assume symmetric nearest neighbour hopping 
of particles with rate u{p, Xi) which implies the following exchange dynamics for {xi}: in time interval dt a length /i 
is transferred from Xi to x^+i with probability u(fj,, Xi)dt/2 and from Xi to Xi-i with probability u{fi, Xi)dt/2. The 
transition is accepted if the constraints ([3]) are obeyed by the updated variables Xi. 

Similarly, we define the nearest neighbour symmetric volume exchange dynamics with rate u'(A, Vi): in time interval 
dt a volume A is transferred from Vi to Vi^i with probability w{A,Vi)dt/2 and from Vi to with probability 
w{A,v^)dt/2. 
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B. Factorised Steady State 

In order to obtain a solvable case that may allow us to study possible condensation scenarios, we seek a factorised 
steady state p3 | where the probabilities of a microscopic configuration {xi, Vi} are of the simple form 



P{xi,vi, . ..xnjVn) = A 



N 



Y[gix^,V,)e{x^ - V^^) 



,i=l 



S{L-Y,^:)S{V-J2v.) (4) 



where g[xi,Vi) are single-particle weights and A is a normalizing constant. That is, the probability of a configuration 
factorises into a product of one factor g{xi, Vi) for each particle. The ^-functions impose the constraint Q at each site 
and the ^-functions impose the global constraints of length conservation and volume conservation. We now assume 
that the single site weight g itself factorises 

g{xi,Vi) = a{xi)b(vi) . (5) 

A sufficient condition for the stationary state to be of the form with g of the from ([5]), is that the x and v dynamics 
independently respect detailed balance with respect to a{x) and b{v). Note that the constraint ([3]) will not enter into 
this requirement of detailed balance since we demand detailed balance between any two of the configurations allowed 
by the constraint. Therefore, we require V/i, x, x' , A, v, v' 

u(ii,x)a{x)a{x') = u{^,x' + ^)a{x — ^)a{x' + ^i) (6) 
w{A,v)b{v)b{v') = w{A,v' + A)b{v^A)b{v' + A) (7) 

which imply that u{fi, x)a{x) / a{x — fi) is independent of x and w{A, v)b{v)/b{v — A) is independent of u, leading to 

u{fi,x) = dfi)^^^ 8 
a(x) 

wiA,v) = d(A)fc^ (9) 
b{v) 

where c(/i) and d(A) are arbitrary positive functions. Therefore rates of the form (|8l9p lead to an equilibrium state 
of the factorised form ^ with single particle weights ([5]). Moreover, appropriate choice of the rates (|8l9p allow any 
single-site weight ([5]) to be generated. 

In the following, for simplicity we restrict ourselves to g{xi,Vi) = 1. This requires — c(/i) and w{A,v) = 

d(A), in order that a(x) = b(v) = 1. In this case the jump size /x of a p-sphere does not depend on the separation 
Xi and A, the amount of volume transferred from one p-sphere, does not depend on the volume Vi of the p-sphere 
(up to the constraint that the remaining volume is non-negative). The condition g{xi,Vi) = 1 implies that all allowed 
microscopic configurations have the same steady state probability. Thus the thermodynamics of the system are driven 
entirely by the constraints in the configuration space coming from the hard-sphere condition i.e. any phase transition 
will be entropy-driven. One can easily generalize our calculation for arbitrary b{v). It turns out that the condensation 
transition occurs for other choices of b{v) as well, as long as the function b{v) decays exponentially or slower for large 
V. However, we stick to the choice b{v) = 1 in this paper for simplicity. 

III. THERMODYNAMICS IN GRAND CANONICAL ENSEMBLE 

A. Ensembles 

The natural ensemble generated by the dynamics discussed in the previous section is the micro-canonical ensemble 
wherein only microscopic configurations with the correct total length and total volume are allowed and all allowed 
configurations have the same statistical weight in the steady state given by 



P{{Xi,Vi}) 



Zn{L,V) 



N 



<5(L-^x,)5(F-^^;0 (10) 



where L and V are given and the normalizing constant Zn{L, V) is just the micro-canonical partition function for N 
particles on a ring given by the integral over the allowed microscopic configurations 

N 
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i.e. Zn{L, V) is a volume in configuration space. We define the two basic control parameters: the density of particles 
p = N/L and the volume per particle (j) = V/N. Our goal in this section is to compute, in the large N limit but 
for given fixed (p, 0), the single-site distribution p{x,v) obtained from the joint distribution (fTO|l by integrating out 
all the {xi,Vi} variables except at one site where they are held fixed with values {x,v). Next, from this single-site 
distribution p{x, v) we will derive the marginals p{x) — p{x, v) dv and p{v) = p{x, v) dx for given (p, (j)). We 
will show that the marginals p{x) and p{v) exhibit a rich variety of behavior in different regions of the (p, (p) plane, 
including a condensation transition for p > 1 . 

To make progress, we shall first follow the usual route of replacing the hard constraints on the allowed total length 
and total volume in the micro-canonical ensemble by the soft constraints of the grand canonical ensemble where 
the total length and volume are allowed to fluctuate. Thus, for this system we define the grand canonical partition 
function as the double Laplace transform of the micro-canonical partition function Zjv(i, V) 



/ dL dVZN{L,V)e 
Jo Jo 



-LX-Vs 



Then inserting ((TTl) into ([12]) yields 



N 

n 

N 



dXjdVr 



) e 



.-1 "'0 



da;,; 



where the single-particle partition function G{X, s) is defined as 



G{X,s) 



due 



-si; — At; 



l/p 



(12) 

(13) 

(14) 
(15) 

(16) 



Thus in the grand canonical ensemble the joint distribution pop factorises P{{xi,Vi}) — YiiPi^ii^i) with the single- 
particle joint distribution p(x, v) given by 



pix, v) = [G(A, s)]-^ e-^^-'^eix - v^/P)e{v) 



(17) 



where we have made explicit the fact that v > 0. One can integrate out x or v to obtain respectively, the single-particle 
volume and length marginal distributions 



p{x) = [G{\,s)]-'- 



i/p_ 



-\x 



e{x) 



(18) 
(19) 



The constraint 9{x — v^^p) in (|17p makes the x and v variables manifestly coupled. Note however that p{x, v) can be 



'diagonalised' or decoupled if one uses the gap variables gi 
p{x,v) =p{g,v) becomes 



Xi 



,1/P 



(see Fig. [ij instead of the Xi variables. Then 



pi9,v) = [GiX,s)]-' [e-^^%)] 



e-'"-^"''" e(v] 



(20) 



The consequence of this decoupling will be discussed later. 

It remains then to estimate the two Lagrange multipliers (A, s) by enforcing the conservation of total length and 
total volume on an average. The mean lattice length L and mean volume V in the grand canonical ensemble are given 
by 



L = 



dlnZffiX, s) 
dX 



V 



51nZjv(A, s) 
ds 



(21) 



In this ensemble, our original control parameters (p, (p) are thus replaced by the averages 



N 



V 
N 



(22) 
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Using the definitions ((2T|) and the expression of Zjsf (fT5|l gives 

1 91nG'(A,s) 91nG'(A,s) 



(23) 



Thus, for given values of the control parameters (p, (/)), we have to solve the two conditions ((23l) to get the corresponding 
values (A, s) and then use them in (jlSp and (|19p to obtain the marginals within the grand canonical framework. 

For the discussion in Section |V] on the alternative free energy functional route, it turns out to be convenient to 
introduce another dimensionless observable denoting the line coverage 

^ = f (24) 

where Lp = N{v\^^) is the average total length occupied by the p-spheres. In other words, rj is the ratio of {v\^^'), the 
average length of a p-sphere, to 1/p, the average available length. The first average is also easy to compute within 
the grand canonical ensemble using ([T| 



AG(A.., -|mAG(A,.,|. (25, 

The right hand side can be expressed in terms of the density (I23p and we get 

A = -^. (26) 
1 — 77 

The quantity on the right hand side of ((26l) is nothing but the pressure of a polydisperse hard-rod fluid [1^. As 
expected, we find that the Lagrange multiplier A associated to the conservation of total length L coincides with the 
pressure (see section IVB I) 

Substituting A from (|26p in (|20p , wc find that the marginal gap distribution p{g) — Jq°° p{g, v) dv for this polydisperse 
hard-rod fluid is given by 



p{g) = exp 
1-77 



P9 



0{9)- (27) 



We note that in the case of monodisperse hard rods (the Tonks gas, where all rods have the same length) the 
equilibrium gap distribution has exactly the same form as ([27)) [2l| . Here we see that expression ([27]) is more general 
and holds even for the gap distribution of a polydisperse hard-rod fluid. 

B. Reduced volume fraction ip*' and scaling variable u 

To facilitate further analysis it is useful to consider a scaling combination of the two Lagrange multipliers s, A 

(28) 

If we make a change of variable v ~ X^p y the single-particle partition function G(A, s) in (|16|) can be written in the 
scaling form 

where 

H{u) = / dye-"^-^'^^ (30) 

Similarly, the constitutive equations for the density and volume fraction (|23p may be cast in a scaling form in terms 
of the variable u 



Fp{t-.) (31) 
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where the two scahng functions Fp{u) and F^{u) are given by 

W = 1 + "roo'' , (33) 



J^dye-'^y-y' 

J^°^dyye--y-y"^ dlng(u) 
^^^"^ - ^-^dye-"^-^^^" ' d^ ■ ^^^^ 

As a result, we can eliminate A between pip and p2p by introducing the scaled volume per particle, 0*, defined by 

0* = (35) 

which depends only on a single scaled variable u — s/X^: 

^ = -F>r ■ ^ ^ 

We call (j)* the reduced volume fraction. Such a reduction to a single scaling variable (instead of two independent 
variables {p,4')) can be traced back to the fact that in the grand canonical ensemble the joint distribution p{x,v), 
when expressed in terms of the gap gi variables, essentially decouples as in (|20p . In due course we will investigate the 
behavior of the function (p* (u) for different values of p. 



C. Entropy 



As already stated the condensation transition in our model is completely entropy driven, therefore we should first 
define the entropy. In terms of the original micro-canonical partition function (jlip. the entropy per particle is 



N 



(37) 



Since the configuration space volume Ziy{L,V) (|lip may be less than 1, S, as defined above, has no reason to be a 
positive definite variable. The entropy could of course be made positive by dividing Zjv by a suitable small constant 
in the same way that the partition function for an ideal gas contains a factor h~^^ . 

In the grand canonical ensemble, the entropy per particle may be expressed as the Legendre transform 



1 



In 2^Ar(A, s) — A 



d\nZi^{\,s) (91nZjv(A, s) 



ds 



lnG(A, 



A 

h sq 

P 



where in the second line we impose the constraints 

dS_ 



. A consequence of ([55]) is that 



-P 



dp 



^ A 



(38) 
(39) 

(40) 



In particular, we note that the derivative of S with respect to vanishes at s = 0. 

We now consider the entropy as a function of the scaling variable u — s/A^, introduced in Eq. ((28| . For fixed 
density p, one can express the entropy 5 as a single function of u. Substituting ([3T|) . ([32| and ([29| in (|39l) allows us 
to write 



S = -{p+l) lnp-(39+l) \iiFp{u)+liiH{u) + Fp{u) + uF^{u). 



(41) 



p varies, u varies according to (|36p and consequently the entropy S 



For fixed density p, as (j) and hence 
varies. 

In Appendix [Al we show from expression (I4ip that the entropy achieves a maximum when u — 0. With this 
knowledge we will investigate in the next two subsections how the entropy S, for a fixed density p, behaves as a 
function of (p* for p < 1 and for p > 1. The results are shown in Figure [21 We will see that a condensation transition 
occurs only for p > 1 where, when (jj* is increased above a critical value, u sticks to the value u — and the entropy 
sticks to its maximal value. A vestige of the condensation transition can also be seen for p < 1 in the fact that the 
entropy decreases when (f>* exceeds the critical value (and u passes through zero) which refiects that the configuration 
space becomes severely constrained. 
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FIG. 2: Entropy plot of S' = S + {p + 1) log p, where 5* is the grand canonical entropy, as a function of reduced density , for 
different values of parameter p. The three corresponding critical densities are <j>* = 2/{3n) ~ 0.212 for p = 1/2, <j>* = 1/2 for 
p = 1 and (j}* = 4(2/7r)^''^/5 ~ 0.688 for p = 3/2. In the first two cases no transition occurs. For p > 1, the plateau seen for 
(/>* > (/)* is a signature of the condensate formation. Note that for p < 1, the maximum allowed value of (j)* is 1, whereas one 
can have </>* > 1 for p > 1. 



D. The case p < 1 



Recall that, according to our model, for p < 1 the volume of the p-sphere increases sub-linearly with its diameter 
and for p = 1 the volume is equal to the diameter. Let us start by making the following observation. By definition, 
the inverse density is l/p = L/N = Due to the constraint Xi > v^^, it follows that l/p > {J2i'^y^)/^- 

For p < I, one can use Jensen's inequality to write 



l/p 

J_ 

N 



= (42) 



Thus, for p < 1, one must necessarily have (f>* = (jy^^^p < 1. In other words, for a given p, one is physically allowed to 
increase (f) only up to i.e., (j) < p~^ ■ 

Now, let us fix the density p and imagine increasing (j) from to its maximally allowed value p"''. In other words, 
the reduced volume fraction 0* = 4>^/Pp increases from to its maximally allowed value 1. For a given 0*, we then 
have to solve for u = s/A^. 

The case p < 1: Consider first the case when p <1 strictly. The case p — I will be discussed subsequently. For p <l, 
we notice that the integrals in ([33]) and (I34|) are convergent for any u G [~oo, oo\. Their leading asymptotic behaviors 
can be easily deduced. Specifically, one finds that as m ^ oo 

FM -> 1 (43) 

F^{u) ^ -. (44) 
u 

On the other hand, as u — > —oo, to leading order, 

Fp{u) -> (45) 
F^{u) ^ iM^/^i-f). (46) 

As a consequence, the function <j>*{u) in pop has the following asymptotic behavior 

(t)*{u) ^ 1 as u^-rx) (47) 
^ u-'^'^' as u^oo. (48) 



9 



In addition, one can check that (j)*{u) is a monotonically decreasing function of u, achieving its maximally allowed 
value 1 as u — > —oo (see Fig. [3]). Thus, for any given 0*, we can always find a solution u to the equation (j)* — 4>*{u) 
where 4>*{u) is given in ([36|) . Knowing this solution u, one finds subsequently A > from ((3T|) and s from the relation 
s = uX^. This means that the grand canonical framework works over the full allowed range < 0* < 1 and the two 
marginals p{v) and p{x) have always the form in (fT8|) and (fT9)) with A and s determined as above. This shows that 
there is no condensation for p < I. 




-10 -5 5 10 



U 



FIG. 3: The function (j)*{u) vs u for p = 1/2. 



However, a vestige of the condensation transition still remains even for p < 1. To see this, imagine again that 
we increase the value of the control parameter 0* — cjj^^Pp from to 1. As 0* is increased, the solution u to the 
equation decreases monotonically from oo to — oo (see Fig. Note that when the solution u hits 0, the 

corresponding value of (^*(0) can be computed by putting m = in (|36p and carrying out the elementary integrals 
that gives Fp{0) = {p+l) and ^^(0) = r{2p)/T{p). Thus 



no) = 



1 



1+p 



r(2p) 



r(p) 



i/p 



(49) 



Now, let us investigate the entropy S in ([55]) as a function of the control parameter 0*. As is increased 
monotonically from (and consequently u decreases monotonically from +oo), one can check that the entropy S, 
expressed as a function of u as in (|4T|) . increases monotonically as long as u > 0, i.e., (f>* < <j>*{Q) given in ([49]). For 
(j)* > </>*(0) (or u < 0), the entropy decreases with increasing (jj* . Thus, the entropy S has a maximum value at 
(j)* = 0*(O), see Fig. [2l This is also evident from (|40|) where the derivative of the entropy as a function of (p (for fixed 
p) vanishes at s = and hence at m = 0. Thus the value of 4> which the entropy becomes a maximum (for a fixed 
p) can be appropriately denoted by 4>max{p) and is given by 



4>max{p) = 



10) 



1 p 



1 1 r(2p) 

pP {l+p)P T{p) 



(50) 



The corresponding maximum value of the entropy is obtained by putting m = in (|4ip . Using -Fp(O) = (p + 1) and 
i/(0) = T{p + 1) from ^ we get 



Smax =p+l + \nT{p+l)-{p+ 1) ln((p + l)p) 



(51) 



Physically this means that for cj) > (pmaxip) or equivalently in terms of the reduced volume fraction, for <j)* > (/'*(0), 
the configuration space becomes constrained resulting in the reduction of entropy. 

Can one see a reflection of this vestige of a condensation transition directly in the volume distribution p{v) in Eq. 
P^ ? Indeed one does observe a change of behavior of p{v) as the control parameter (/)* increases through (/'max(p). 
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As (j)* increases from 0, the solution u of = u is positive as long as (/)* < (t>max{p) (see Fig. |3]). Consequently 

s — u\P is also positive and hence the distribution p{v) oc exp[— Au^/^ — sv\ is a monotonically decreasing function 
of V with its maximum at t; = 0. However, when 0* > 4'max{p), the solution u and hence s becomes negative. 
As a result, oc exp[— Aw^/'P+|s|i'] now develops a maximum at a nonzero characteristic volume f* — (p|s|/A)P/^^~^-*. 

The case p = 1: The general conclusion reached above for p < 1 remains valid even for the marginal case p = I, though 
the details are slightly different. Indeed, for p = 1, we can obtain explicit solutions for the marginal distributions. 
To see this, let us first explicitly express A and s in terms of p and </). The integrals in p3p and p4p can be easily 
performed explicitly for p — I giving 

W = ;i±|; (52) 

Thus, unlike the p < I case where the allowed range of u was u G [—00,00], for p = I, u has the allowed range 
M G [—1,00]. Eq. (l36|) then gives 

riu)^^. (53) 
u + 2 

Thus (j)* {u) is again a monotonically decreasing function of u in u G [— 1 , 00] and for any given (jj* one can always find 
a solution u — l/4>* — 2. Hence, as in the case p < 1, there is no condensation for p = 1 as well. 
Note that at u = 0, 0* (0) = 1/2 and hence (j>max{p) = l/(2p). Also, (pij) and ([5^ yield explicitly 

^ (54) 



s + A 

- - - Y ' (55) 
P A 

determining A = p/(l — p(t>) and s = (1 — 2p0)/ (1 — p(/)) in terms of p and 0. This also gives, using (|24p . rj = p(f) — (f>* . 
The grand canonical joint distribution for the separation and volume of a particle (jl7p is 

p(a;,i;) = A(s + A)e-^^-^"6l(j:-i;) . (56) 

The marginal distributions (|18ll9p become 

p{v) = ^ (57) 

P(^) = , ^ , [e-P"/(^-P^) -e-"/^j . (58) 
1 — 2 p L J 

Thus, although there is no phase transition for p — 1, there is a change in the leading behaviour oip{x) as the volume 
per particle </> is increased past <f>max{p) = l/(2p): for < l/(2p) the large x behaviour is p{x) ~ e~^''/^"'^~'"'^^p/(l — 
2p(j)); for (j) — l/(2p) the large x behaviour is p{x) = xe~'^^'^/(jP (note that this expression, for (j) — l/(2p), holds 
for all values of left-to-left distances x); for > l/(2p) the large x behaviour is p{x) ~ e~^/'^p/(2 0p — 1). Thus for 
> l/(2p), the exponential decay of p{x) is the same as that of p{v). 

The entropy per particle ([55)) reads 

5 = ln(l-p0)-lnp + ln0 + 2. (59) 

Note that, as mentioned below the entropy S is not restricted to be only positive. This is evident in the p = 1 
case from (I59p where S can be negative for certain values of the parameters p and 0. One can verify easily that the 
entropy S in (j59p has a maximum at = 4>max{p) = l/(2p) with value 5max — 2 + ln(2) — 21n(2p) from (j5ip . see 
also Fig. [21 As increases past l/(2p) the entropy decreases, implying a constrained configuration space. Thus the 
scenario for p = 1 case is similar to p < 1 discussed earlier, as seen clearly in Fig. 2 for the two representative cases 
p = 1/2 and p = 1. 



E. The case p>l 

For p > 1, the integrals on the right hand side of ([55)1 and ([M)) are convergent only for w > 0. Thus the lowest 
allowed value of u is 0. When u ^ 0, the function 0*(w) in ([5S)) approaches 0*(O) which is still given exactly by (I49p . 
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Thus this is the maximum value of 4>* allowed by the grand canonical ensemble. If 4>* exceeds </'*(0), u = s/X^ cannot 
decrease below 0. It then sticks to its value u — and all the extra volume condenses into a single site. Thus (/)*(0) 
is the critical value beyond which the grand canonical ensemble breaks down signalling the onset of a condensation 
transition. 

At this critical point, the volume fraction = [(f)*(0)/p\P will again be denoted by (pmaxip) and has the same 
expression as in the p < 1 case, namely 

. N 1 1 r(2p) , , 

pP {l+p)P T{p) 

Also, letting m = in (PT|) gives the critical value = {p + l)p. Consequently, the two marginals in and at 
the critical point become 

p{v) ^ [G{0,X)]-'\e-'^''' = j;^^ [p(l+p)f e-^^i+rt^^^^ (6i) 

p{x) ^ [G{0,X)]-'e~^^xP^—^[p{l+p)]'+PxPe-P^'+P^^ (62) 

Thus at condensation p(v) changes from (dominant) exponential decay (|18p to a slower stretched exponential decay 
([M]) and p{x) changes from an exponential decay ([TO]) to the exponential decay multiplied by x^ ([S^ . We interpret 
these results as describing the critical fluid; in the condensed phase we expect a condensate to coexist with the critical 
fluid. For later purposes, we also note that at the critical point the dimensionless line coverage 

V^Vc = (63) 
[l+p) 

which follows by substituting Ac = (p + l)p in (|^. 

How does the entropy S behave as (j) increases from to 4'max{p) for fixed p? As (f) increases monotonically, u 
decreases monotonically till u hits 0. Consequently, the entropy S in (|4ip increases monotonically up to u = 0, 
achieving a maximum value at u = (or equivalently at = (j^maxip)) given by the same expression as in the case 
p < 1 in (j5ip namely, 

S^ax=P+l + lnT{p+l)-ip+l) H{p+l)p)). (64) 

What happens to the entropy when </) exceeds the critical value (/'max(p)j i-e., when a condensate sets in? To see 
this, we note that 4>max{p) in (|60p can be neatly expressed in terms of the p-ih moment {x^) of the critical marginal 
p{x) in ((62|) . This moment can be easily computed and comparing to the expression of (pmaxip) in (|60p one easily 
verifies that 

Ct>max{p) = l(xn ■ (65) 

The physical meaning of (|65p is that condensation occurs when the volume per particle is equal to half the mean 
available volume per particle in the grand canonical distribution. In the case p < 1 this is the value of above which 
the entropy decreases. For p > 1, the intuitive explanation is that at this point, rather than the entropy decreasing 
as was the case for p = 1, the entropy can be held constant by one particle containing the excess volume and leaving 
the rest of the system at the critical volume fraction, see Fig. [51 

Let us then summarize the p > 1 case. The existence of the grand canonical solution means that the system is in 
a fluid state. Thus, for fixed p, 0max(p) hi (|60p is precisely the critical line in the (p, 0) phase diagram, as shown in 
Fig. (ID) for p ^ 2. When the volume fraction cj) exceeds this critical value 0max(p), the grand canonical description 
breaks down and a condensate forms in the system. The entropy S increases monotonically with increasing (j) till (j) 
hits the critical value (pmaxip) where it achieves its maximum value Smax given in (j64p . When 4> is increased to a 
value (j) > 4>max if) , the entropy remains at the maximal value through one particle containing the excess volume and 
forming a condensate. 



F. The large p limit 

The limit p ^ oo is worth mentioning also as it makes links with some other well studied problems. As mentioned in 
the begining of Section 2, in the p oo limit our polydisperse hard-rod system reduces to the classical monodisperse 
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FIG. 4: Phase diagram in the {p,(f>) plane for p = 2. The critical line 4>max{p) = 3^ (for p — 2) separates the fluid and the 
condensed phase. 



Tonks gas of hard rods where each rod has the same size 1. In this case, our result for the gap distribution in pT)) 
(valid for general p), reduces to that of the classical Tonks gas. 

On the other hand (again for p — > 00), the volume of each rod becomes a "passive" property (a color v € [0, oo[). 
From (jlSp . this property is exponentially distributed {p{v) oc exp(— sw)). Such a Poissonian distribution may have 
been anticipated, and is readily obtained within the free energy functional formalism presented in section |Vl We note 
that this exactly coincides with the exponential distribution of wealth and income obtained in the model presented 
in reference [2^ . 

We also note that in the large p limit, the condensation transition disappears. Thus the system is always in a fluid 
state and the volume statistics obey "ideal gas" behaviour. In the hard core language of section |Vl we will indeed see 
that the critical line fraction tends to unity for large p, so that the parameter range (/) > 6rnax actually corresponds 
to an unphysical region where "hard rods" necessarily overlap. Hence, the model of Ref [24] shows no condensation 
transition (more complex interactions between the "agents" would be required). 



IV. MICRO-CANONICAL ANALYSIS OF THE CONDENSED PHASE FOR p > 1 



For p > \ the grand canonical ensemble can only support a volume per particle < (f>maxip), therefore in order 
to fully analyse the condensed phase, where (p = V/N > (pmaxip), one needs to work within the micro-canonical 
ensemble. 

To compute the micro-canonical partition function one inverts the Laplace transforms in (|15p 



ci+00 



C2+00 



dA 



sV + \L 



(66) 



where ci and C2 are chosen so that the integration contours are to the right of any singularities. Using the expression 
rsl) one finds 



Zn{L,V) 



C2+00 



dA 
27ri* 



(67) 



where 



ipi\,s) =^ - + sq 
P 



In A + In 



(68) 
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The integral (|67|) may be evaluated by the saddle-point method. The saddle-point equations coming from the condi- 



tions II = and If = read 



1 dv v^/Pe-'"-^" 

A ^ dv e— 

dv e-^-^-''" 



i/p 



(69) 
(70) 



which are, of course, precisely the grand canonical equations pil [32)1 . Thus when the saddle point exists (i.e. in the 
fluid phase) the results of the grand canonical and micro-canonical ensembles coincide. However, in the condensed 
phase 4> > 4>max{p) one can no longer solve the saddle-point equations for s > 0. Therefore Ziq{L^V) must be 
evaluated by an alternative approach. 

It will be useful to consider the Laplace transform of Zn{L, V) with respect to the length L (rather than the double 
Laplace transform of Zn{L, V) which generates the grand canonical partition function). 



(71) 
(72) 



Zn{X,V) = / dLe-^''ZN{L,V) 
Jo 

^ />oo />OC 



N 



■ N 

n 



dvj{vt 



S{J2v,-V) 



(73) 



where we have defined 



(74) 



This definition ensures that the integral of f{vi) is normalised to unity therefore f{vi) may be considered as the 
probability distribution for a positive random variable Vi. Then the quantity 



f^jv(A) 



■ N 

n 

.i=l 



dvtf{vt) 



(75) 



manifest in (j73p . is the probability that a sum of N independent positive random variables each distributed according 
to / is equal to V. Assuming that / has finite first and second moments /ii and /i2, as it does in the case (j74p . where 



P-r, 



dv v''f{v) 



(76) 



one can invoke some limiting results on the sum of a large number N of such random variables. We first define 

Vc = Nn, (77) 
then the following results for sums of random variables derived in a different context jT^] will be useful 

1 (V-Vc)^ 



For V-Vc^ 0{N^/^) f7jv(A, V) 
For V >Vc and V ~ Vc 0{N) 



nN{X,V)^Nf{V^Vc) 



(78) 
(79) 



where A'^ — ^2 ~ Pi- The first result is a central limit theorem which expresses the fact that the sum is Gaussian 
distributed about the mean Vc- The second is a large deviation result whose interpretation is that for the sum of 
random variables to be equal to a value V, much greater than the mean Vc, one of the random variables should be 
equal to V — Vc to leading order and the other — 1 should be of 0(/ii). In ([75]) the factor f{V — Vc) comes from 
the probability of the large random variable and the factor N comes from the number of ways of choosing the large 
contribution from the N random variables. 
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Let us note that fii is precisely the critical volume fraction 4>max{p) defined in the previous section. This follows 
by computing the first moment of f{v) in (j74p and comparing it to the expression of (jimaxip) in (|60p . Thus 



(80) 



justifying the subscript c (critical) for the volume V . 

We may now obtain forms for Zm{L^V) near to criticality and in the condensed phase by inverting the Laplace 
transform 



Ci +2C30 



Z.iL^V) = I ^e 



LA 



A 



N 



(81) 



Here the A integral may be evaluated without problem using the saddle-point method with the saddle located (to 
leading order for large L, large N but keeping the ratio p = N/L fixed) at A* = (p + l)p. For rijv(A*, V) one can use 
the results in ^ and This gives 



Zn{L,V) 



Zn{L,V) 



N 



2Tr{p+l)N 



2'k{p+ 1 

1/2 



1/2 



pr(p) 



pr(p) 

(A*)P+i 



e 2JVA'= 



Af-1 



(A,)P+i 



for F - Vc - 0(iV) 
for V-Vc^O{N^/^) 



Therefore for > (f>rnax{p) the entropy per particle, given in the micro-canonical ensemble by 



remains fixed at 



S=-\nZN{L,V) , 



s = p + 1 + \iiT{p + 1) - {p + i)Iti Hp + i)p) 



(82) 



(83) 



(84) 



(85) 



Note that this is precisely the maximal value Smax in (|64p computed in the fiuid state. To summarize, for p > 1, 
the entropy S increases monotonically with increasing (f), achieves its maximal value Smax at = </imax(p) a-nd then 
remains fixed at this value for all > (j)max{p)- 

We may also consider the marginal distribution p{v) in the condensed phase. As a signature of condensation this 
distribution should contain a bump at v = V — Vc- In the micro-canonical ensemble we have 



p{x,v) 



Zn-i{L ~x,V -v) 
Zn{L,V) 



i/p) 



(86) 



Integrating out x yields 



p{v) = Zn{L, Vy^ / da; Zn-i{L - x,V - v) 



i/p 



(87) 



One can then substitute the asymptotic behavior of the function Zn{L, V) obtained in ([82ll and (|83ll to estimate p{v) 
in ([57]) . Omitting details, we find that 



For V <t:V ~Vc p{v) 



A* 



r(i+p) 



For V = V -Vc + 0{N'^'^) p{v) 



1 



1 



N V27r7VA2 



("-(V-V"c))^ 
= e 2NA'^ 



(89) 



The latter piece (15^1) of ^(w) represents the condensate: it has a total weight 1/A^ signifying a single condensate site 
and shows a Gaussian distribution of the bump around the excess volume V — Vc- 
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V. THE FREE ENERGY FUNCTIONAL ROUTE 

Having established that our problem, defined in terms of dynamical rules, admits a factorized steady state prob- 
ability with detailed balance, we can envisage the steady state as the equilibrium state of a hard-rod model, where 
the distribution of rod lengths I = v^Ip is not known o priori. This size distribution should be that which minimizes 
the free energy of the system (or equivalently maximizes the total entropy, since we only deal here with excluded 
volume interactions). Our goal in the remainder is therefore to provide a different, perhaps more physical -but fully 
equivalent" perspective, and obtain the optimal polydispersity minimizing the free energy functional of a hard-rod 
system, with the same global constraints as in section |lll fixed density p — N / L, and fixed p-moment of the length 
distribution {{v) = {t^)), where p is a parameter of the model, and is fixed as in previous sections. In doing so, we 
are only interested here in the size distribution (variable w), and not in the joint distribution of size and gap x. In 
this respect, the approach of sections IlIII and Hvl provides more detailed information. 

The problem we face is the one-dimensional analogue of the optimal packing of polydisperse hard sphere fiuids 
addressed in Refs. pTi . [l9l . [25j . which exhibits an unexpected condensation transition. The free energy functional 
approach put forward here is akin to that used in [T9| . with nevertheless the interesting feature that in the present ID 
geometry, exact results can be obtained. On the other hand, approximations of Percus-Yevick type or more involved 
treatments were necessary in Refs [13, [H, HI] . Indeed, the problem boils down to finding the size dependence of the 
chemical potential of a given species in a polydisperse mixture, which is not known exactly for hard discs or hard 
spheres. 

A. Free energy functional for polydisperse hard rods 

We work in the canonical ensemble where a distribution of N hard rods at temperature T occupies an available 
length L. Their total length is denoted Lp and defines the line coverage (or packing fraction) 77 = Lp/L = p{tj where 
p = N/L is the density and {£) is the mean size. To establish the connection with the p-spheres discussed earlier, we 
can view each rod of size £ as a (hyper)sphere constrained to move on a line, and having "volume" v = . 

Writing the ideal entropy of a multicomponent and discrete system is straightforward (26j . but the limit of a 
continuous distribution requires some care (see the discussion below). For hard rods, the free energy functional may 
be written as a sum of ideal and excess contributions [27j : 

(}T{W} Jd£ W{e) [In {A^pWie)Mi)) - 1] + PT,^{W}. (90) 

where A is an irrelevant length scale, /3 is the inverse temperature, and W is the length probability distribution 
function (such that / di £W{£) = 60)- The excess free energy for a homogeneous system of hard rods is known 
exactly, /3jFox{W^} = — ln(l — 77) [2g|, and can be generalized to inhomogeneous situations [l^. In Eq.f [TO)) , the 
function je{£) -yet to be specified- ensures that different choices for the labelling of the particles will lead to the same 
optimal length distribution W{£). One could indeed choose v — £p (or say any other power) as a working variable, with 
an associated "labelling" function jv{v) and a probability distribution function such that W{£) — p £p~^Wv{£p)- 
Enforcing the consistency of both descriptions imposes ji{£) = £P~'^jv{£^) (up to an irrelevant prefactor), where the 
factor £P~^ is the Jacobian of the transformation £ £p. The natural labelling for the particles, i.e. that which gives a 
constant function j, follows from the way polydispersity is sampled (see e.g. [l^ [lE|)j and reflects the dynamics of the 
system. In the following Monte Carlo simulations, we will attempt to change the size of two particles £1 and £2 selected 
at random by adding a small increment to £1 and conversely subtracting the same quantity from £2, in accordance 
with the rule specified in section [Til and more precisely with the requirement that b{v) = 1. This corresponds to a fiat 
function j^, and hence to = £p~^. However, given an arbitrary b{v) (not constant) in the microscopic rate in Eq. 
([5]) , it is far from evident how to choose an appropriate labelling function j„ (v) in the macroscopic description in Eq. 
(|90)) . that would precisely correspond to this microscopic rate. 

Minimization of J- — with the two constraints of normalization J d£W{£) — 1 and fixed p-moment J d££PW{£), that 
can be included by two Lagrange multipliers — leads to the functional form of the optimal size distribution W* , which 
should then coincide with the probability distribution function of v-^^p, where the statistics of v is given by Eq. (|18p . 
Since the explicit expression of J- is known, see above, this provides a first angle of attack to our problem. However, 
interesting information follows from more microscopic considerations and scaling arguments, as becomes clear below. 
Note that the constraint that the line coverage should not exceed unity is not explicitly considered, but is implicitly 
encoded in the excess contribution to the free energy functional ([50)1 . 
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B. Optimal polydispersity distribution 

We take advantage of the fact that any infinitesimal change in W* has a vanishing free energy cost 6J^, provided 
the two aforementioned constraints are satisfied. We proceed in two steps: a) a given rod of arbitrary length £q is 
expanded : io ^ + ^^o! b) all particles are rescaled {i ^ at) in order to fulfil the constraint of a conserved p-th 
moment. In the first step, this moment changes by an amount 



^ ^ ^p^r''54 (91) 

while in the second stage, it changes by {d^ — 1)(^^) — p{oi — 1){P'). Enforcing the conservation of {P) then imposes 



" - N{IP) > 
In stages a)+h)^ the distribution function changes by an amount 

5W=^ [5(^-4-54) - <5(f-M] + (l-«)^^- (93) 

It proves convenient to express the variation of the ideal part of the free energy in terms of the probability distribution 
function Wv of v. This yields, with vq = £q 

5T;d = kbT < — — — - + — \ 5vQ. (94) 

The excess contribution variation -which is most conveniently expressed in terms of € rather than v- is in addition 
exactly the reversible work required to perform the changes under consideration. In step a), we note that this work 
is the same as if a confining "wall" would be displaced a distance 5^^^ thereby compressing the system. We therefore 
have 

6Ta ^ PSeo, (95) 

where P is the pressure of the system. This relation is s pec ific to the one dimensional case, and at the root of 
important simplifications as compared to higher dimensions [23|. On the other hand, it is a general result that in any 
space dimension, the work performed during step b) is (see e.g. appendix B of [l9i] ) 

STb - -—5L„ (96) 
V 

where Pcx = P ^ pkT is the excess pressure and SLp is the total volume change of the particle upon the rescaling 
l-^ at. SLp^ {a - l)Lp = {a - l)N{e). Making use of Eq. we finally have : 

Combining this with Eq. ((94|) supplemented by the requirement that = for = W* (or equivalently = W*) 
gives 



W*ivo) (v) ?7 p{iP) 



= 0. (98) 



After integration with respect to Wq, we arrive at 

(99) 



K(.)«exp{-/3P.V.+ (^^_l) 



This is, expectedly, the very same form as obtained in section Hill see Eq. (|T8)) . In addition, Eq. (|99)) above makes 
explicit the connection between the chemical potentials A and s appearing in Eq. psp and intensive thermodynamic 
quantities, for example A = (3P, see Eq. ([^S)) . Equivalently, in terms of the £ variable, we obtain 

W%e) = AF~^ exp{~^P£ + - l) i!^} (100) 
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where ^ is a normalization factor. As alluded to earlier, it is noteworthy here that the pressure is exactly related to 
the density p through [H [H, [H 

PP = where r} = p{t). (101) 

1 — 77 

One can note that the low density behaviour of W* allows for some consistency test of our prediction. When the 
density vanishes, one has P and pPcx/v 0, so that 

W*i£)^e-'e^p(^-^y (102) 

hence an exponential distribution oi v = (the conserved quantity). Such an expression also immediately follows 
from direct minimization of the functional ([501) . restricted to its ideal contribution (jF^x — 0). We recover here the 
expression obtained in section Hill in the large p limit. As mentioned earlier, we also recover the results reported in 
an econophysics context 24] for a simple model where agents act as our ideal particles, exchanging random amounts 
of a quantity v (money) in binary encounters. 

With a weight function such as (|102p . one immediately finds 

as it should (the quantity (P) appearing on the l.h.s. of p02p is therefore indeed the moment of order p of the 
distribution). This shows the consistency of our distribution function in the low density limit. The corresponding 
mean size follows, assuming again the low density form (|102p : 

f^£W*(£)d£ 

= r(l + l/p) {£P)^/P, (104) 

With the standard "functional route" alluded to earlier (i.e. direct minimization of Eg. ([50]) ) which does not provide 
explicitly Lagrange multipliers, the requirement (|103p together with normalization would determine those multipliers 
for any density. 

In the remainder, the star superscript will be omitted to refer to the optimal distribution W*, without ambiguity. 



C. The condensation transition 



We revisit here the condensation transition brought to the fore in section 3, in the more liquid-state language of 
hard rods. Since the excess pressure following from (jlOip fulfils the identity Pex/v = P, as can be readily checked, 
the distribution pOOp can be rewritten more explicitly as 



W{£) = A£P~^ exp 



hi"' + 




£P ' 


1 - ?7 




W). 



(105) 



with rjc — p/{l+p). For p < 1, this distribution is normalizable for all line coverages rj (given that 77 < 1). This is no 
longer the case for p > 1, provided that 77 > rjc, where 77^ will be referred to as the critical line coverage. For p > \ 
and r] > r]c, the divergent behaviour of distribution (jlOSp at large £ is indicative of the formation of a macroscopic 
aggregate with size Lq. The scenario is identical to that inferred from the observations of Ref. O and worked out 



in [17| . The system relaxes by transferring "volume" to the aggregate, which effectively acts as a piston and coexists 
with a polydisperse mixture Ai confined in a region of size L — Lq. In this region, the size distribution obeys Eq. 
()105p . where 77 is no longer the total line coverage fraction, but should be replaced by the line coverage 77^ in the 
region free of aggregate. The different line coverages in the problem are connected through 



1 - ?/o 

where 770 — Lq/L is the line coverage of the aggregate. We show in Appendix IB] that 

for rj < rjc 
for rj > r]c- 




(106) 



(107) 
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The line coverage rjo may be considered as an order parameter for the transition. The key ingredient to arrive at Eq. 
P07p is that the fluid outside the condensate (assuming the latter species forms) is critical, in the available length 
L — Lq. It therefore has line coverage ijm = r/c, and size distribution 



W{i) cx £P-'^ exp 



1 - ??c 



(108) 



We note finally that the case p oo is specific in that we then have rjc ^ 1 (while cf)* — > 4/e). Since the hard 
core constraint imposes ?7 < 1, we see here why the large p limit reduces in fact to the ideal gas case, where £p is 
exponentially distributed and no transition occurs. 

VI. COMPARISON WITH SIMULATION RESULTS 

In order to verify analytical predictions it is useful to compare with numerical results. For example, numerical studies 
of condensation have in the past revealed important information about when the asymptotic behaviour predicted 
analytically actually emerges in a finite system. 



A. Control parameter and critical point 

To analyze in more detail the scenario at work and put our predictions to the test, we perform Monte Carlo 
simulations. We first have to introduce a relevant control parameter (we emphasize that except when p = 1, the line 
coverage or packing fraction rj (j24p is not a conserved quantity, but is self-consistently determined). Since both the 
density and {£p) are conserved variables, we will use the dimensionless reduced volume fraction 



(109) 



introduced in cq. (I35|) . A simple convexity argument shows that for p > 1, > 77, while the reverse holds for p < 1. 
Having chosen {£p)^/p as our relevant length scale, we also introduce a reduced pressure 



PP{£P) 



'p\i/p _ 



and a rescaled length 



which, from pOS]) . has size distribution 

W{I) = A'\ 



exp 



l-ry l-p{£)' 

I^£/{£py/p 



_n_J_ ^ f y-Vc 



{£y 



(110) 



(111) 



(112) 



It is in general not possible to relate explicitly the reduced density (/)* to the line coverage 77, except at the critical 
point rj ~ rjc, since then W takes a pure exponential shape (up to the algebraic prefactor). From (jlOSp . we have there 



so that 



Starting from 



v=nc 



Tjp + k) 
(2^-1)!^^/" 



for fc > 



(p-iy. 

when ?7 = 0, the reduced density 0* increases with rj and reaches the value 



1 /(2p-l) 



1+p V (p-1)! 



i/p 



(113) 



(114) 



(115) 



when r] = Tjc- This corresponds exactly to the threshold obtained in section HVl see Eq. ([SO)) . For p > 1, i]c signals 
the onset of the formation of the condensate. For p < 1, rj > rjc only signals a region where the sub-dominant term in 
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P in the exponential (|105p changes sign. The size distribution may therefore change from the unimodal shape found 
at low line coverage to a bimodal form when rj exceeds some threshold, itself larger than rj^. On the other hand, the 
onset of bimodality for the distribution Wy oi v = P is rjc- 

When no condensate forms (i.e. p < 1 or cj) < (jjc H p > 1), we compute numerically the probability distribution 
function of i as follows. For a given value of the reduced density (j)*, the mean value {£) appearing on the rhs in 
()112|1 is determined self consistently by enforcing that it should coincide with the first moment of the distribution 
having statistical weight (|112p . Note that (/>* = rj/{l). Alternatively, we may also impose the second self-consistency 
requirement that {P) = 1. We have systematically checked that both routes provide the same result for from 
which the different quantities of interest may be computed. This illustrates the consistency of the functional form 
(|112l) . With such a procedure, the solution found numerically is always unique; it will be compared against the results 
of Monte Carlo simulations in the next section. 

In situations where a condensate is expected (</> > <j)c and p > 1), an explicit prediction between the control 
parameter (jf and the condensate size (or more precisely condensate line coverage) can be derived from the remark 
that the fluid phase outside the condensate is critical. In other words, 

" \ iGfluid / 

On the other hand, the global reduced density (including, thus, the condensate) reads 

V iGfiuid / 

The two above equations allow us to compute rjQ, the condensate line coverage, through 

0*^^ = iVf-^g + [0:(l-r;o)f . (118) 

Hence, for large N, r/o —* 0, and we have 

^'o - (119) 



B. Monte Carlo simulations 



To test our predictions, we have implemented Monte Carlo simulations, closely following the algorithm used in Refs 
p^ . N hard rods with different sizes are confined on a line of length L; a first type of move amounts to randomly 
selecting a particle, and randomly translating it. A second kind of move allows the system to relax its size distribution, 
and sample polydispersity. Two particles are selected at random in the system; the size £i of particle 1 is expanded at 
the expense of the size of particle 2, so that + £2 is constant: — > £f + A; £^^ £^ — A, where the increment A is 
drawn from a distribution w such that the typical value is small compared to . Both types of moves are accepted 
provided they do not lead to any overlap between the rods and for the second kind, provided the shrunk rod does not 
have a negative length. These rules correspond to the model defined in section [TTl in particular to the case where the 
volume exchange rate ^(A, v) of section 2 does not depend on v. 

Typical results are shown in Figs. [5] and [HI that are both for cases without condensation. In all figures, the data 
gathered from the Monte Carlo simulations arc shown with the symbols, while the predictions are shown by the 
continuous curves. The agreement theory /simulations is very good (see in particular the distribution function in the 
upper inset of Fig. O. From Eq. (|104p . we have that in the low density limit, {£) = 2 for p = 1/2 (see the lower inset 
of Figure (5]), and {£) — \/tt/1 — 0.886 for p = 2 (which can indeed be seen in the inset of Figure 0. 

We now turn to the cases where a condensate should form. It may then be difficult, from a practical point of view, 
to distinguish such a big particle from others belonging to the tail of the size distribution. Another difficulty comes 
from the fact that the condensate line fraction 770 niay be small for large N . Attention must be paid to the size of 
the condensate that is expected. Eq. (|119p indicates that the condensate line coverage scales with N hke A^^i+^/P; 
however, the typical rod size in the fluid phase concomitantly exhibits a faster decay in 1/A'^, so that there is always 
a clear separation condensate/fluid. In our simulations, we have followed for p > \ the biggest particle (size Lq) in 
the simulation box (line). As can be seen in Fig. [71 the corresponding line fraction 770 = Lq/L closely follows our 
prediction pisp . We therefore conclude that Fig. [7] proves the existence of the condensate for (j>* >(/>*. 
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FIG. 5: Equation of state as a function of line coverage 77 in the case where the exponent defining the conserved moment is 
p = 1/2. The symbols are for the Monte Carlo results, and the curve is for Eq. (|110|l . A characteristic length distribution 
is shown in the upper inset for r) ~ 0.72 (corresponding to 0* = 0.6). The distribution (|112|) is compared to its Monte 
Carlo counterpart. The lower inset shows the first moment {^) as a function of line coverage. Here, rjc — 1/3 for which 
02 = 2/(37r) ~ 0.212 and {£)c = 7r/2. Monte Carlo data with A*' = 1000 particles are shown with the circles. 
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FIG. 6: Equation of state as a function of line coverage 77 when p = 2, in which case the critical line coverage rjc — 2/3 corresponds 
to (j)* — y^2/3 ~ 0.816. The exact pressure is shown by the continuous line, while the dotted curve is for the crude approximation 
where {£} is assumed density independent and equal to its critical value -^2/3, which leads to P* — -\/3/2 r]/{l — rj). Circles : 
Monte Carlo results with A'' — 1000. The inset shows the density dependence of {£). 

To confirm the results of section 4 where the condensate volume distribution is calculated, see Eq. ([55]) . we measured 
the fluctuations in the condensate volume. We find that that the distribution p{v) turns out to be Gaussian, as 
predicted in section 4 (see Fig. [5]). Another prediction is that the fluid phase outside the condensate should have a 
size distribution of the form (|108[) . irrespective of the reduced volume fraction provided 0* ></)*. This expectation 
is fully consistent with the Monte Carlo results, see Fig. [9] for p — 3/2. For convenience, we have considered in the 
main graph the distribution of £p, since it is predicted to be a pure exponential. We have performed a similar analysis 
at p = 5/2, for various volume fractions beyond the critical one, and the same conclusion holds with again a critical 
fluid phase, and a size distribution (excluding the condensate) that does not depend on the imposed volume fraction 
(j). The fluid phase can only accommodate a well deflned flnite fraction of the total volume (or length), so that when 
(/)* increases above the extra volume is transferred to the aggregate (see Fig. [7]). 
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FIG. 7: Packing fraction (line coverage) of the condensate as a function of reduced volume fractions (j>* , for p = 3/2 (lower 
sets) and p = 5/2 (upper sets). The symbols have been obtained by computing the mean size of the largest particle observed in 
Monte Carlo simulations, and the curve shows Eq. (|118p . For the present parameters (A'^ = 1000), the approximation provided 
by Eq. (|119p proves quite accurate (dashed line, hardly distinguishable from the continuous curve for p = 5/2). 
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FIG. 8: Linear-log plot of the aggregate size distribution for p — 3/2 and two reduced volume fractions (0* = 0.9 and 2). The 
data are obtained following over time the different sizes taken by the (unique) condensate. Here, P denotes the probability 
distribution function of y = v'' /{Vo) ~ 1- 



VII. CONCLUSION 



We have considered a simple stochastic model of mass transport where a real-space condensation takes place upon 
increasing the volume fraction. The system studied consists of p-spheres, constrained to move on a ring, with hard 
core constraints : the left to left particle distance Xi between p-spheres i and i + 1 should exceed the diameter £i of 
p-sphere i (see Fig. [T]). The latter quantity is itself related to the volume of the sphere through ii — v'J^ , where 
p has been retained as a parameter. The total volume of the particles, 'Ylii'^i-> is fixed, so that the present model 
involves two conserved quantities (^ ■ Xi, the total length available, and the total volume V = Vi). The dynamics 
comprises diffusion (which can be seen as a stochastic exchange of the quantity x between neighbouring particles), 
and a stochastic exchange of volume. When p > 1 and beyond a critical density that has been worked out explicitly, a 
particle of macroscopically large "mass" (large volume) appears, carrying a finite fraction of the system total volume 
surrounded by a critical fluid phase, with a size distribution independent of the total density. 

Under mild conditions pertaining to the sampling of the conserved quantities, we have shown that the system 
admits a factorized steady state probability density, with detailed balance between allowed configurations (essentially 
those with no overlaps between the p-spheres) . This allows us to use an alternative description of the system in terms 
of finding the equilibrium distribution of an ensemble of hard rods, that do not have a quenched size distribution, but 
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FIG. 9: Linear-log plot of the probability distribution function Wv of particle "volumes" v = £^ as a function of £, for p — 3/2 
and several reduced volume fractions beyond the condensation transition (circles for (j)* — 2, squares for (j)* = 5 and diamonds 
for <;/>* = 8). Here 0* = 4(2/^)^/^5 ^ 0.688. The continuous curve shows the prediction of Eq. (|108|l . The inset shows the 
probability distribution of the variable £ instead of £^, again as a function of ^ = v^^^.) 



can exchange length provided the global constraint £^ remains fixed. Formulated as such, the problem appears as 
the one dimensional analogue of the optimal polydispersity studies of hard discs and hard sphere fluids [13, EM HE] 
(for instance, the problem studied in [191 co rresponds to the three dimensional hard sphere situation with p = 1). 
The phase transition reported in [itI . llSl. Il9l [25| can therefore be seen as a condensation arising from constraints in 
the configuration space. 

We have characterized the scenario at work for this transition, and obtained analytically the probability distributions 
of volumes, gap distance g, and left to left particle distance x. Restricting to the volume distribution, we have shown 
that a density functional approach supplemented with scaling considerations allows to recover the results derived from 
the stochastic processes viewpoint (working out the consequences of the factorization property of the steady state). 
In this respect, the problem is to find the polydispersity of a hard rod system that will minimize its free energy, given 
that the total number of rods is fixed as that the p moment of the diameter distribution is fixed as well. The consistent 
predictions of both approaches have finally been successfully tested against Monte Carlo simulations. It is remarkable 
that this optimality problem leads to a phase transition in a ID system, where the condensate size can be considered 
as an order parameter (note that the constraints imposed introduce a global coupling between all particles). It should 
also be emphasized that whereas the higher dimensional systems in two or three dimensions could exhibit a transition 
for p — 1 (with the convention that the conserved quantity for hard spheres in d dimensions is (^''^)), such a value 
for p-spheres on a line turns out to be critical with however no transition observed: a macroscopic aggregate can only 
form provided p > 1 . 

An important question that remains is that of the dynamic pathways to condensation, that is, how does a single 
condensate emerge from some given initial condition for the p-sphere volumes. As the condensation is driven by 
constraints in configuration space the dynamics too might be strongly affected by constraints, possibly producing, for 
example, entropic barriers. This issue remains to be explored. 
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APPENDIX A 

In this appendix we show that the grand canonical entropy S{u) defined in (|4ip as 

S{u) = -{p+1) \np-{p+l) \n Fp{u)+ In H{u) + Fp{u)+uF^{u). (Al) 

is maximised when the scaling variable u takes the values zero. 

Expanding about w = to second order in w = and using the definitions of H{u), Fp{u) and F^{u) given in 
section 3.2 we obtain 

S{u) = (1 + p - (1 + p) ln(l +p)+ In [F(l + p)]) + am + a2U^ + 0{u^) (A2) 
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where 

ai = (A3) 
02 = [{l+p + p'')T\2p)^{l+p)T{p)T{3p)]/[2{l+p)T\p)] (A4) 

The first equation shows that u = is a stationary point. The function 

cip) - (1 + p + p2)r2 (2p) - (1 + p)Tip)Ti3p) (A5) 

is maximised at p = 0.614413 where its value c{p) — —0.566304 is negative. Hence c{p) is negative for all p > 0, which 
implies 02 < for all p > proving indeed that at u = 0, S{u) is a maximum. 



APPENDIX B 



To find out what does determine 770 (and thus rjj^ ) for a given value of 77, we reconsider the free energy functional 
([SO]) , supplemented with two Lagrange terms to fulfil the constraints [13) EM HE] 

TZiW) = (3T{W) + NCo I W{e)d£ + NCp ( [ PW{1) Al + LU . (Bl) 



We have to minimize this expression with respect to W and Vb- Stationarity with respect to W leads to an expression 
of the form (fT05)l : 

\og[l^~m{l)] + + £0 + C,P = 0, (B2) 

while the derivative with respect to Vq reads 

dn d[3T p_i 

The free energy of the system is a function of L — Lq, since a constituted aggregate does not contribute to F apart 
from the confinement it induces on the remaining particles which have free energy : F{N, Lq, L) = Fj^ {N, L — Lq). 
Hence, Eq. (|B3p also reads 



P + NpCpLP-\ (B4) 



From pOSp and (jB2p . we have £p cx {rjc — rjj^). On the other hand, a physically acceptable mixture M., should have a 
normalizable size distribution, which imposes 77^ < 77c- Hence, £p > so that the derivative (|B3|) cannot vanish, and 
is always positive (a situation already encountered in [l3|)- The minimization with respect to Lq therefore leads one 
to choose for Lq the minimum possible value compatible with 77^^ < rjc- The optimal value of is thus for rj < r]c, 
and such that 7;^ — rjc whenever rj > rjc- In other words 

{0 for r] < rjc 
^ for 77 > 77.. (B5) 
1 - 7?c 
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